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NON-UNIFORM  MAGNETOHYDROEVNAMIC  SHOCK  PROPAGATION, 
WITH  SPECIAL  REFERENCE  TO  CYLINDRICAL  AND  SPHERICAL  SHOCK  WAVES 


Roy  M.  Gundersen 

Summary 

Previous  work  on  non-uniform  shock  propagation  in  monatomic 
conducting  gases  is  generalized  to  an  arbitrary  value  of  the  adiabatic 
index.  Specifically,  the  perturbation  generated  when  an  initially  uniform 
hydromagnetic  shock  of  arbitrary  strength  impinges  on  an  area  variation  is 
determined,  the  problem  being  linearized  on  the  basis  of  small  area  variations. 
When  the  shock  encounters  the  area  change,  the  shock  strength  is  altered,  and 
the  subsequent  flow  is  non-isentropic.  There  are  two  distinct  contributions  to 
the  perturbation,  namely,  a  permanent  perturbation  due  to  the  area  change  and 
a  transient  reflected  disturbance,  and  expressions  for  these  are  obtained. 

A  first  order  relation  between  area  change  and  shock  strength  is 
obtained  and  integrated  numerically  to  give  an  area- shock  strength  relationship 
valid  for  channels  with  finite  continuous  area  variation.  Particular  area 
distributions  are  utilized  to  discuss  converging  cylindrical  and  spherical 
hydromagnetic  shocks. 

The  present  theory  includes  results  on  non-uniform  gas  dynamic  shock 
propagation  as  a  special  case. 

Sponsored  by  the  Mathematics  Research  Center,  U.  S.  Army,  Madison,  Wisconsin, 
under  Contract  No.  DA-11-0E2-ORD-2059. 
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1.  Introduction 

In  recent  years,  there  has  been  a  great  deal  of  work  on  shock  propagation 
in  channels  with  small  non-uniform  cross-sectional  area  distributions.  By  the 
use  of  a  linearization  based  on  small  area  variations,  Chester  [  1]  found  that 
the  pressure  perturbation  produced  when  an  initially  plane  shock  wave  passed 
through  a  non-uniform  transition  section,  which  joined  two  channels  of  constant 
but  unequal  cross-sectional  area,  was  given  by; 

-KiP^  -  Pj)[AA]/A, 

where  P_  -  P,  was  the  Initial  pressure  discontinuity  across  the  shock,  [A A] 
the  net  change  in  area  and  the  parameter  K  a  monotonlcally  decreasing  function 
of  the  shock  strength  with  only  a  small  total  variation. 

Chester  started  with  the  full  three-dimensional  equations  of  motion  and 
then  carried  out  an  averaging  process  by  considering  only  the  average  pressure. 
This  final  restriction  means,  however,  that  the  same  results  must  be  obtainable 
by  investigating  the  problem  from  the  very  beginning  by  a  one -dimensional 
approach.  This  very  important  simplification,  which  permits  of  great  extension, 
was  first  observed  by  Paul  Germain,  and  Chester's  results  obtained  through  a 
one -dimensional  analysis  by  Gundersen  [6],  extracts  from  which  appeared  in 
[7  ],  the  analysis  being  based  on  characteristic  perturbations  presented  by 
Germain  and  Gundersen  [5]. 

It  is  surprising  that  the  simple  one- dimensional  approach  led  to  an 
improvement,  namely,  the  term  [  AA]/A  in  Chester’s  work  was  replaced  by 
^A  where  A  was  the  perturbed  area  distribution.  This  showed  clearly  the 
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first  order  differential  relation  between  the  area  change  and  shock  strength 
at  any  point  in  the  non-uniform  section.  The  significance  of  such  a  relation 
was  first  realized  by  Chisnell  [  2] ,  who  was  led  to  it  in  a  different  way, 
essentially  by  assuming  that  Chester's  steady  state  solution,  valid  for 
large  time  only,  could  be  utilized.  Chisnell  used  an  integrated  form  ( in 
closed-form)  of  the  shock  strength-area  relationship  to  give  an  approximate 
description  of  the  motion  of  the  shock  in  terms  of  the  area  of  the  duct.  By 
suitable  choices  of  the  area  distribution,  a  description  of  converging  cylindrical 
and  spherical  shocks  was  given,  and  the  results  checked  by  comparison  with 
previous  similarity  solutions,  valid  in  the  neighborhood  of  the  points  of 
collapse  of  the  shocks.  The  comparison  showed  the  remarkable  accuracy  of 
Chisnell's  work. 

Extensions  of  Chisnell's  work,  including  closed-form  shock  strength- 
area  relations,  for  piston-driven  shocks,  where  reflected  waves  from  the 
piston  come  back  to  interact  with  the  shock  and  modify  its  strength,  were 
given  for  linear  and  quadratic  area  distributions,  the  appropriate  ones  for 
cylindrical  and  spherical  shocks,  in  [8]  and  [9].  Those  results  are,  of 
course,  applicable  to  reflections  from  a  contact  layer. 

Piston-driven  shocks  propagating  through  ducts  with  arbitrary,  though 
small,  area  variations  were  considered  by  Mirels  [15]  by  a  source  distribution 
method  and  in  [10]  by  the  small  perturbation  method,  and  the  solution  given 
in  infinite  series  form. 

The  problem  originally  considered  by  Chester  has  been  extended  to 
the  case  where  the  fluid^n  front  of  the  shock  was  not  at  rest  in  [11]  and 
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to  the  case  that  the  area  perturi^atlon  was  also  time 'dependent  in  [  12] . 

Recently,  many  others  have  worked  on  similar  or  related  problems, 
e.  g. ,  in  [16]  and  [17],  Chester's  results  were  rederived.  Recently,  it 
was  shown  [13]  that  in  a  monatomic  conducting  gas  subjected  to  a  transverse 
magnetic  field,  the  perturbation  produced  when  an  initially  uniform  plane 
hydromagnetlc  shock  of  arbitrary  strength  encountered  an  area  variation 
could  be  determined  the  problem  being  linearized  on  the  basis  of  small 
area  variations,  and  the  solution  presented  in  a  form  which  included  the 
usual  gas  dynamic  results,  described  in  the  foregoing,  as  a  special  case. 

In  (14),  a  theory  rather  parallel  to  that  of  Chisnell  [2]  was  developed  and 
converging  cylindrical  and  spherical  hydromagnetlc  shocks  discussed.  The 
limitation  of  a  monatomic  gas  was  due  to  the  fact  that  analogs  of  the  usual 
Riemann  Invariants  could  then  be  determined  explicitly.  In  the  present 
paper,  it  is  shown  that  this  restriction  is  unnecessary  in  spite  of  the  fact 
that  analogs  of  the  Riemann  invariants  cannot  be  determined  explicitly,  at 
least  in  terms  of  elementary  functions. 

The  problem  considered  is  that  of  an  initially  uniform  plane  hydromagnetlc 
shock  travelling  with  constant  speed  into  an  ideal  gas  at  rest  in  a  duct  which 
has  a  section  of  non-uniform  area.  To  ensure  that  all  changes  in  the  motion 
of  the  shock  are  due  to  the  area  variations,  it  is  assumed  that  these  are 
confined  to  the  region  x  >  0,  whereas  to  the  left  of  the  cross-section 
X  =  0,  the  tube  is  of  constant  cross-sectional  area,  and  the  shock  propagates 
with  initially  constant  speed  in  this  portion  of  the  tube. 

When  the  shock  encounters  the  area  variation,  it  is  perturbed,  the 
shock  strength  altered  and  the  subsequent  flow  non-isentropic.  There  are 
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two  distinct  contributions  to  the  disturbance,  viz. ,  a  permanent  perturtiation 
due  directly  to  the  area  variation  and  a  transient  disturbance,  due  to  reflections 
from  the  shock  of  the  permanent  disturbance,  which  propagates  with  velocity 

r  2 . .ZA  • 

CO  =  [c  +  b  , 

where  c  is  the  local  speed  of  sound  and  b  the  AlfvSn  speed,  with  respect 
to  the  flow  behind  the  shock. 

In  the  neighborhood  of  the  shock,  the  ultimate  effect  is  an  altered 
shock  strength  and  concomitant  pressure  change  behind  the  shock.  When 
the  main  flow  behind  the  shock  has  speed  <  co ,  the  transient  reflected 
disturbance  is  convected  to  the  left  with  speed  w  relative  to  the  fluid 
while  for  flow  speed  >  co ,  this  disturbance  is  convected  to  the  right 
with  speed  co  relative  to  the  main  flow.  This  latter  must  be  added  to  the 
steady  flow  solution.  Expressions  for  these  various  contributions  are 
obtained. 

Specifically,  it  is  shown  that  the  pressure  perturbation  immediately 
behind  the  incident  shock  is  given  by: 

where  perturbations  are  denoted  by  a  bar,  is  the  original  pressure 

discontinuity  across  the  shock,  the  cross-sectional  area,  a  the  shock 
strength  defined  as  a  density  ratio  and  m^  =  ,  a  measure  of  the  applied 

transverse  field. 

For  all  m  ,  lim  K  =  0.  5  and  11m  K  =  0.  39414  for  y  =  7/5  . 

(r-l+  0—6-  ^ 
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But  these  are  precisely  the  limits  for  the  corresponding  parameter  in  the  gas 
dynamic  case  which  corresponds  to  m^  =  0  so  there  is  the  result:  For  very 
weak  or  very  strong  shocks  [  in  an  asymptotic  sense] ,  the  results  are  independent 
of  the  applied  field  and  agree  with  the  usual  gas  dynamic  results. 

In  magnetohydrodynamics,  strong  shocks  occur  for  <r  close  to 
(-Y  +  1)/(Y-  1)  ,  where  y  Is  the  adiabatic  index,  or  for  a  very  strong 
applied  field  for  any  o-  >  1  ,  i.  e. ,  is  then  large. 

For  =  0 ,  is  a  monotonically  decreasing  function  of  the  shock 
strength  and  agrees  exactly  with  the  corresponding  parameter  in  the  afore¬ 
mentioned  publications.  For  any  o,  there  no  longer  is  a  monotonic 

variation  with  <r  ,  but  each  curve  is  concave  upward  with  curves  for  greater 
m^  lying  beneath  those  for  lesser  and  all  curves  pass  through  the  points 

(o-,  K^)  =  (1>  0*5)  and  (6,  0.39414)  for  y  =  7/5  .  Further,  for  fixed 
incident  shock  strength,  is  a  monotonically  decreasing  function  of 
m^  ,  1.  e. ,  in  a  diverging  (converging)  channel,  the  pressure  decrement 
(increment)  is  decreased  (decreased)  by  increasing  the  applied  field. 

If  the  area  variations  are  confined  to  a  transition  section  of  finite 
extent  joining  two  portions  of  constant  cross-sectional  area,  the  total  effect 
of  the  passage  through  the  transition  section  is  an  altered  shock  strength 
with  the  shock  asymptotically  becoming  again  uniform.  Qualitatively,  the 
motion  of  the  shock  is  independent  of  whether  the  main  flow  behind  the 
shock  has  speed  <1^  or  >  . 

Finally,  the  first  order  differential  relation  between  shock  strength 
and  area  variation  is  integrated  numerically,  and  particular  area  distributions 
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are  used  to  consider  converging  cylindrical  and  spherical  magnetohydro- 
dynatnic  shocks  which,  in  contrast  to  plane  shocks,  are  inherently  unstable 
and  ultimately  become  strong.  This  section  generalizes  the  results  of  [  14j 
to  arbitrary  values  of  the  adiabatic  index  y,  though  numerical  tables  are 
presented  only  for  y  =  V5 . 

Specifically,  it  is  shown  that  the  strengths  of  converging  cylindrical 

and  spherical  magnetohydrodynamlc  shocks  near  the  centers  of  these  fronts 

-K  -2K 

are  proportional  to  D  and  D  ,  respectively,  independent  of  the 
applied  field,  where  D  is  the  distance  from  the  centers  and  K  =  0.  39414 
for  y  =  7/5 . 

The  theory  of  gas  dynamic  shock  propagation  in  non-uniform  ducts 
Is  contained  as  a  special  case  of  the  theory  presented  in  this  paper,  i.  e. , 
it  is  obtained  simply  by  setting  m^  =  0 ,  so  that  there  is  a  check  on  the 
theory  presented  herein. 
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2.  General  Theory 

The  quasi-one-dimenslonal  non- steady  motion  of  an  ideal,  inviscid, 
perfectly  conducting  compressible  fluid  subjected  to  a  transverse  magnetic 
field,  1. e. ,  the  induction  B  =  (0,  0,  B),  is  governed  by  the  system  of  equations: 

2c^/{  Y  -  1)  +  2ucy  { Y  -  1)  +  +  ucA^A  =  0  (  2. 1) 


u^  +  uu  + 

t  X 


2ccy(Y-l)  +  b^ByB-c^s^/Y(Y-l)c^ 


=  0 


(2.2) 


B,  +  uB  +  Bu  =  0 

t  X  X 


(2.3) 


s^  +  us  =  0  (2.4) 

t  X 

2  2  - 

where  u,  c,  s,  p,  b  =  B  /|xp,  p.,  y  ^nd  A  are,  respectively,  the  particle 
velocity,  local  speed  of  sound,  specific  entropy,  density,  square  of  the 
Alfv^n  speed,  permeability,  ratio  of  specific  heat  at  constant  pressure  c^ 
and  at  constant  volume  c  and  cross-sectional  area  of  the  channel. 

V 

Partial  derivatives  are  denoted  by  subscripts,  and  all  dependent  variables 
are  function  of  x  and  t  alone  save  A  which  is  considered  time  independent. 
The  characteristics  of  this  system  are: 


where 


+  b 


—  =u,  u  +  w,  u-W 


:-r 


V 
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For  an  arbitrary  isentroplc  constant  area  flow,  the  base  flow  in  the 
neighborhood  of  which  perturbations  will  be  considered,  the  basic  system 
of  equations  { 2. 1) ,  ( 2.  2)  and  (  2.  3)  may  be  written  in  the  following 
characteristic  form  [  3] : 

Xp  -  (u  +  w)tp  =0  (2.  5) 

X  -  (u  -  w)t  =  0  ( 2.  6) 


-  ut^  =0  (  2,  7) 

-  c^)Bp/B  +  wup  +  2ccp/(y  -  1)  =  0  (2.8) 

-  c^)B  /3  -  wu  +  2cc  /(y  -  1)  =  0  (2.9) 

a  a  a 


B^/B  -  2c^/(Y  -  l)c  =  0  (2.10) 

with  characteristic  parameters  (a,  p,  |). 

From  equation  (2.10),  it  is  clear  that  the  quantity  B/c^^^  "  is 
constant  along  each  particle  path,  which  is  a  well-known  consequence  of 
the  assumption  of  Infinite  electrical  conductivity,  i.  e. ,  the  magnetic  field 
is  "frozen"  into  the  fluid,  and  for  a  constant  state  or  a  simple  wave  flow, 

B/c^^^  is  constant  throughout  the  flow  so  that  equations  (2.8)  and 


(  2.  9)  may  be  written  as; 
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i 

i 

s 

f 

i 


Up  +  2wCp/(  Y  -  l)c  =  0 

(2.11) 

-u  +  2«c  /(y  "  l)c  =  0 

a  a 

(2. 12) 

Since  p  =  B  =  with  constants  r^  and 

b^=r^V/<^-Vr3  =  kc^/<^-') 

so  that 

=  c^  [1  +  kc^] 

where 

\  =  2(2-y)/(Y-1)  . 

For  the  case  of  a  monatomic  gas  y  =  5/5  ,  and  equations  (  2. 11) 

and 

(2.12) 

may  be  integrated  explicitly  to  yield; 

u/2  +  [  1  +  kc]  =  u/2  +  ( w/c)  Vk  =  a 

(2. 13) 

-  u/2  +  [l  +  kc]^^Vk  =  -  u/2  +  (w/c)  Vk  =  P 

(2.14) 

where  (a,  P)  may  be  considered  as  generalizations  of  the  usual  Rlemann 
invariants. 

For  an  arbitrary  value  of  y,  the  solutions  of  equations  (10)  and  (11) 
may  formally  be  written  as  ; 


r  X  - 

u/2  +  J[l  +  kc  dc/  (  Y  -  1)  =  a  (  2. 15) 

-  u/2  +  /  [  1  +  kc^]  *  dc/ (  Y  *  1)  =  P 


(2.16) 
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The  solution  of  the  present  problem  for  a  monatomic  gas  was  obtained 

by  a  transformation  of  dependent  variables  which  made  ( 2. 13)  and  (2. 14) 

linear  relations.  This  same  technique  is  employed  for  an  arbitrary  value  of 

y.  Although  the  integrals  appearing  in  equations  (2.15)  and  (2.16)  may 

be  evaluated  in  terms  of  hypergeometric  functions,  an  explicit  evaluation  is 

r  X  i 

unnecessary.  The  introduction  of  the  new  dependent  variable  w  =  j  [1  +  kc  ]*dc 
transforms  equations  (2.1)  and  (2.2)  to: 

w^/(  Y  -  1)  +* uw^(  Y  -  1)  +  «u^/2  +  uwA ^2A  =  0  (  2. 17 ) 

u^/2  +  uu^2  +  ww^(Y-I)  =  c^sy2Y(Y- l)c^  (2.18) 

=  2c ^(y-1)c  and  (2.3)  omitted. 
Adding  and  subtracting  equations  (2.17)  and  (2.18)  gives  the  system  to  be 
utilized  in  the  sequel: 

w^/(Y-1) +u^/2+(u  +  w)[u^2  +  w^(Y-1)]  =-uwA^2A+c^s  ^2y(y-1)c^  (2.19) 
w^/(  Y-1)  -  u^/2  + (u  -  w)[ -u^/2  + w^(  Y-1)  ]  =-uwA^2A-c^s j^2y(Y-1)c^  (2.20) 


where  use  has  been  made  of  the  relation  B 


s^  +  us  =  0 

t  X 


(2.  21) 


The  close  analogy  of  this  system  with  the  basic  equations  of  gas  dynamics 
should  be  noted,  viz. ,  a  hyperbolic  system  of  three  non-linear  first  order  partial 
differential  equations  with  Riemann  invariants  which  are  linear  relations  between 


-12- 


#310 


the  dependent  variables.  The  system  is,  in  fact,  the  magnetohydrodynamic 
generalization  of  such  and  includes  the  gas  dynamic  equations  as  a  special 
case.  For  k^O,  use,  w=c  and  the  system  (2.19)  -(2.21)  reduces 
exactly  to  that  used  by  Gundersen  [  6],  [  7]  to  obtain  the  solution  of  the 
propagation  of  a  shock  in  a  tube  of  varying  cross-section  and  perturbations 
of  simple  waves.  The  above  system  has  also  been  utilized  to  obtain  the 
perturbation  of  a  magnetohydrodynamlc  simple  wave,  and  the  results  will 
be  presented  in  a  forthcoming  paper. 

For  the  present  problem,  it  is  only  necessary  to  perturb  in  the 
neighborhood  of  a  constant  state  (isentropic),  and  a  formal  linearization  of 
equations  (2. 19)  -  (2.  21)  in  the  neighborhood  of  this  known  state,  denoted 
by  the  subscript  zero,  leads  to  the  following  system  of  linear  equations  for 
the  terms  of  first  order,  denoted  by  the  subscript  one: 


=lt  + 


u-s  =  0 
0  lx 


(2.  24) 


where  R  =  u^/2  +  w^/(-y-1)  ,  S  =  -u^/2  +  w^/('y-1)  . 

According  to  (2.24),  s^  remains  constant  along  the  particle  paths  of 
the  given  flow,  i.  e. ,  along  dx/dtsu^.  Since  Pq(  dx  -  u^dt)  is  the  exact 
differential  of  a  function  4^^  which,  when  equated  to  a  constant  defines  the 
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particle  paths,  the  solution  of  (2.  24)  may  be  written  as: 

Sj  =n(4'o) 

with  n  an  arbitrary  differentiable  function. 

It  is  convenient  to  define  a  new  function  Tp(x,  t)  by 

Vlx  =  Vo“''V ■ 

Then,  the  general  solution  to  {2.  22)  -  (2.  24)  may  be  written  as; 

R  =  F[x- (uQ  +  WQ)t]  +E[x-UQt]/wQ  -  UqW^Aj/2Aq(Uq  +  «q)  (2.25) 

S  =  G[x- (uQ-w^)t]  +  E[x-UQt]/a)Q  -  UqWqAj/2Aq(Uq  -  w^)  (2.26) 

Tq  =  2E‘[x  -  u^t]  (2.27) 

in  terms  of  three  arbitrary  functions  of  one  argument.  From  the  general  solution, 
the  four  distinct  contributions  to  the  perturbation  may  be  noted,  namely,  one 
due  to  the  entropy  variations,  which  travels  along  the  particle  paths  and  is 
measured  by  E ,  a  disturbance  due  directly  to  the  area  variations,  a  perturbation 
propagating  with  velocity  ,  the  true  speed  of  sound,  with  respect  to  the 
fluid  along  the  family  of  characteristics  x  -  (^Q  +  w^lt  =  constant  and  measured 
by  F  and  a  perturbation  propagating  with  velocity  with  respect  to  the  fluid 

along  the  family  of  characteristics  x  -  (u^  "“q)!  =  constant  and  measured  by  G. 
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3.  Solution  in  the  Vicinity  of  the  Incident  Shock 

Jump  conditions  across  normal  hydromagnetic  shocks  have  been 
considered  by  several  authors,  e.  g. ,  Friedrichs  [4] .  Let  U  be  the  shock 
velocity,  v  =  U  -  u  and  the  subscripts  one  and  two  designate  flow  quantities 
in  the  regions  in  front  of  and  behind  the  shock.  Then  the  analogs  of  the 
Ranklne-Hugoniot  relations  are: 

Pl'^l  =  P2^2 

BiVi  = 

Pl^'l  "  P2''2  +  ^2  + 

YP/{Y-1)P^  +  Vj/2  +  =  yP2/(Y-1)P2 

For  gas -dynamical  shocks,  knowledge  of  the  flow  in  front  of  and  one 
parameter  behind  the  shock  suffices  to  give  a  complete  solution.  For  magneto¬ 
hydrodynamic  shocks,  all  parameters  behind  the  shock  may  be  expressed  in 
terms  of  those  in  front  and  two  parameters,  e.  g. ,  the  shock  strength  and  one 
which  gives  a  measure  of  the  applied  field. 

Let  <r  =  92/9^  >  ^  >  n  =  u/c ,  M  =  v/c  and 

2  2  2  2 

q  =w/c  =l  +  m  .  Then  the  following  relations  hold,  where  n^  is  assumed 


zero  and 


0  =  (Y+1)/(Y-1): 
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Vj/v^  =  p/Pi  -  B/Bj  =  =  T 

T  =  [Otr-l  -  vin^(l-(r)^/2]/(e-(r) 

c/Cj  =  tA 

m^/nij  =o-Vt  (3.1) 

=  {mj[Y/(Y-l)  +  (2-Y)<r/(Y-l)]  +  2/(  y-l)  }<r/{  O-cr) 

=  M^A  T 

=  (l-o-’^M^lorA)^ 

which  express  the  flow  parameters  behind  the  shock  in  terms  of  <r  and  m^  . 

The  effect  of  an  area  variation  on  the  motion  of  an  initially  uniform 
hydromagnetic  shock  propagating  with  constant  speed  into  a  fluid  at  rest 
will  now  be  determihed  by  the  use  of  the  general  solution  for  the  non-isentropic 
perturbation  of  a  constant  state  as  given  by  equations  (2.25)  -  (2.  27)  .  The 
particular  method  of  generation  of  the  shock  is  left  open  save  that  what  ever 
that  may  be,  the  mechanism  is  sufficiently  far  removed  so  that  no  reflections 
come  back  to  interfere  with  the  basic  interaction  considered  herein.  From 
the  formulation  of  the  problem,  there  is  no  mechanism  downstream  of  the  shock 
which  could  give  rise  to  the  term  F[x  -  ( )t]  in  (2.25)  so  that  the 
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pressure  perturtsatlon  behind  the  shock  may  be  obtained  by  setting  F  =  0 . 

This  gives: 

u/2  +  w^/Cy-l)  -  (3.  2) 

where  perturbations  of  a  base  quantity  are  denoted  by  a  bar.  Since 

^2  =  Y-l)c^  =  V^^^Z  ■  ’ 

the  latter  from  the  equation  of  state,  equation  (3.  2)  may  be  written  as: 

;/2  +  +  '=2V^“2l''’2  ‘  '  ‘'2“2V“2'“2  *  “2' 

or 

^Z^'^z  [  Y(m2  +  ^)‘^JP2^V(Y-l)P2<i2  "  "  ""z^Z^Z^^  "z  ‘^Z^^Z 

(3.3) 

From  the  jump  conditions,  equations  (3.1) : 

(0^-l)(r  +  m^^Yo-[3e-l-60(r  +  3(e+l)(r^- 2<r^J/2  I-  • 

(0-o-)[e(7-l- Ym^(l-(r)  V2] 


/o-' 


2(0+1)  e/Y+  2(0+l)(  0-2)(r/Y  +  m^[  0(l+0)+(l+0)(  0-2)<r+{  0-3)(  20-l)(r^+(  3-0)0- 
2(  0-O-)  (o--!)  {2{  0+l)/Y  +  m^[  (1+0)  +  o-(0-3)J} 


SO  that  equation  (3.  3)  gives; 
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^2/^2  “  " 


(3.4) 


where 


K,  = 


(e^-l)(r  +  m^Yo-l 30-1-  600-  +  3(0  +  l)<r^  -  2<r^]/2 
00-  -  1  -  ■Ym^(l-<r)V2 


m2(0-<r) 

q2(v-i)  ^ 


2(0+1)  0/V+2{0+l)(0-2)(r/Y  +  mj^{0(l+0)+(l+0)(  0-2)0-  + (0-3)(20-l)<r^+(3-0)<r^} 
2{(r-l){2(0  +  l)/Y  +  m^^[{l+0)  +o-(0-3)]} 


*  ^ 

Ylm^+D-l 

(0^-l)o-  +  mjYo-{30- 1-600-  +  3(  0  + l)o-^-2o-^}/2 

v(Y-i)q2 

00-  -  1  -  (l-o-)^/2 

-1 


or 


P2  =  -K2(.,mj)(P2-Pj)A/A2 


(3.5) 


where 


•^2  = 


Ym^(l-(r) V2  +  1  -  0<r 
(0+l)(l-(r)  +  Ymj{l-(r)V2 


1 


Strong  shocks  can  occur  in  hydromagnetics  in  two  ways,  namely,  for 
<r  close  to  (y-^ll/lY-l)  >  e.  g. ,  6  for  y  =  7/5  ,  or  for  a  very  strong 

applied  field  for  any  o-  >  1 ,  i.  e. ,  m^  is  then  large. 

For  all  m  and  y  ■  7/5 ,  lim  K  =  0. 5  and  Urn  K  =  0.39414;  but 

(r-l+  (r-6- 

these  are  precisely  the  limits  for  the  corresponding  parameter  in  the  gas  dynamic 
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case  which  corresponds  to  =  0  so  there  is  the  result: 

For  very  weak  or  very  strong  shocks  [  in  an  asymptotic  sense] ,  the  results 
are  independent  of  the  applied  field  and  agree  with  the  usual  gas  dynamic 


results. 


For  =  0 ,  is  a  monotonic  ally  decreasing  function  of  the  shock 
strength  and  agrees  precisely  with  the  corresponding  parameter  in  the  gas 
dynamic  theory.  For  any  m^^  0  ,  the  monotonicity  is  lost,  but  each  curve 
is  concave  upward  with  curves  for  greater  m^  lying  below  those  for  lesser 
mj^  ,  and  all  curves  pass  through  the  points  (cr,  K^)  =  (1,  0.  5)  and 
(6,  0.39414)  for  y=V5.  Further,  for  fixed  incident  shock  strength, 
is  a  monotonically  decreasing  function  of  m^  so  that  in  a  diverging  (converging) 
channel,  the  pressure  decrement  (increment)  is  decreased  (decreased)  by 
increasing  the  applied  field.  Qualitively,  the  motion  of  the  shock  is  independent 
of  whether  the  main  flow  speed  behind  the  shock  is  <  w  or  >  w  . 

Graphs  of  are  presented  in  Figures  1 ,  2,  3  and  4  for  y  =  VS . 
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4.  The  Reflected  Disturbance 

To  complete  the  solution,  an  expression  for  the  reflected  disturbance 
must  be  obtained  and  this  is  done  through  equation  ( 2.  26).  The  evaluation  of 
the  arbitrary  function  G  is  most  readily  accomplished  by  noting  that  the 
system  (2.  25)  -  (2.  27)  may  be  written  as: 

u^/c^  +  =  -  n n^  +  q2) 


2- 


Thus  on  addition,  which  serves  to  eliminate  u  : 

c» 


2  2 
(V-D  +  Yrn^  m^ 


(e-v)[eo-  -l-Ym^^(l-v)V2] 


q2A2(Y-l)  |(e^-l)v  +  m^ Yo-[3e-l-6ev  +  3(0  + 1) v^-2(r ^ J/2 


2  - 
"2‘^2^2 
/  2  2^  ^ 

2  '  *^2^^2 


=  G[x- (u^-w^ltj/c^ 


Evaluating  this  on  the  shock,  x  =  U  t ,  and  replacing  P^/P,  by  its 

w  Ct 

value  from  (3.4)  gives: 


G[k]/C2  =K^A2[6V]/A2 


where 
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\  bx-(u2-«2)^ 


K  s 

2  2 
*^2”^2 


6  s  (M2+  *^2^/^^2  ■*■  **2^ 


'1 


{Y-i)q, 


(Y-1)  +Ym2  ni2(e-<r)[  e«r-l-  Vm^il-o’lVz] 


{0^-l)(r  +  m^Y<^[3e-l  -60(r  +  3(0+ l)(r^  -  2or^]/2 


Thus: 


~  ^1^2^ ^^“^^2  '  ^2^2^  ^^^'^^2 


(4.1) 


where; 


e,  = 


"2*^2 
2  2 
V'^2 


0<r-l  -  Yni^( l-<r )V2 
(0  +  l)((r-l)  +  Ym^^(cr-l)V2 


1  1^2 
Y-1+ Ym2 


m. 


(0-0-)  [  00- -1  -  Ynij  (l-(r)V2] 


q2Y(Y-l)  (Y-l)q2  j(e^-l)(r  +  injY‘>-[30-l-60(r  +  3(0  +  l)(r^  -  2(r^]/2 


e2  =  K2-^ 


The  parameter  6  is  a  monotonic  increasing  function  of  cr  (for  all  m^)  and 
varies  between  the  limits 

0.  5  <  6  <  1.  6429 

for  all  mj  and  y  =  7/5.  Graphs  are  presented  in  Figure  5.  It  might  be  noted 
that  Chester's  graph  of  the  corresponding  parameter  is  incorrect,  and,  apparently, 
the  reciprocal  of  the  true  value  has  been  plotted. 
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When  the  flow  speed  behind  the  shock  is  <  ,  a  disturbance  will 

be  reflected  to  the  left  downstream  of  the  shock,  and,  from  (4.1),  the 
pressure  perturbation  will  ultimately  be  given  by: 

since  =  0  for  x  <  0  . 

For  flow  speed  >  behind  the  shock,  the  reflected  disturbance 
will  be  convected  to  the  right  and  the  pressure  perturbation  is  given  by  (4.1), 
Graphs  of  the  parameter  ^re  presented  in  Figure  6  for  y  =  7/5. 
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5.  The  Shock  Strength  Area  Relation 

Equation  (3.  5)  is  a  first  order  differential  relation  between  area  and 
shock  strength,  namely, 


dA/A+ dr/iT-DK^Co-,  m^)  =  0  (5.1) 

where 

T=[0o--1-  Ym^(l-cr)V2]/(e-<r)  , 
and  the  solution  may  be  written  as: 


A 


V  exp 


[-/ 


dr 


(T-l)K^(a-  ,  m^) 


(5.2) 


or 


A  =  V  F((r  ,  m^) 


(5.3) 


where 


log^  F((r  ,  m^)  =  -  / 


0^-1  +  m^  y[5  6-1  -  60(r  +  3(l  +  0)a^-2(r^]/2 


(  0-(r)(o--l)  [  0  +  1  +  Ym^((r-l)V2]  K2((r,  m^) 


do- 


and  V  is  a  constant  of  integration.  The  integral  is  given  in  terms  of  cr  since 
is  known  in  terms  of  o-  . 

Equation  (5,  Z)  contains  the  results  of  Chisnell  as  a  special  case,  viz. , 
m^  =  0  and  may  be  considered  as  the  magnetohydrodynamic  generalization  of 
such.  Chisnell  showed  that  an  integrated  form  of  the  shock  strength-area 
relation  could  be  utilized  to  discuss  channels  with  finite  continuous  area 
variations  and  obtained  a  check  on  his  theory  by  choosing  particular  area 
distributions  and  showed  that  converging  cylindrical  and  spherical  shocks 
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could  be  discussed  with  remarkable  accuracy.  Equation  (5.  2)  also  includes 
as  a  special  case  previous  work  on  monatomic  conducting  fluids  presented  in 
[14]. 

Some  comments  can  be  made  immediately  about  converging  cylindrical 

and  spherical  hydromagnetic  shocks  which  can  be  discussed  as  particular 

cases  by  choosing  linear  and  quadratic  area  distributions.  The  cylindrical 

shock  corresponds  to  a  wedge-shaped  channel  with  cross-sectional  area 

proportional  to  the  distance  from  the  center  of  the  shock  front  while  the 

spherical  shock  corresponds  to  a  conical  channel  with  area  proportional  to 

the  square  of  the  distance  from  the  center.  Although  plane  shocks  are  stable, 

converging  cylindrical  and  spherical  shocks  are  unstable  and  as  the  shock 

converges,  its  strength  Increases  and  ultimately  becomes  singular  at  the 

center.  In  the  neighborhood  of  the  singular  point,  (5.  2)  shows  that  A 

,  ,  * 

*  1/ K  ^ 

is  proportional  to  t  ^  ,  where  K  is  the  asymptotic  limit  of  ,  so 

that  the  shock  strengths  of  converging  cylindrical  and  spherical  hydromagnetic 

*  # 

-K  -2K 

shocks  near  the  center  are  proportional,  respectively,  to  D  and  D  , 
where  D  is  the  distance  of  the  shock  from  its  axis  or  point  of  symmetry. 

But  K  =  0.394141,  independent  of  the  applied  field,  and  this  is  exactly 
the  value  Chisnell  obtained  for  gas  dynamic  shocks.'  Thus  the  result: 

Near  the  axis  or  point  of  symmetry,  the  strengths  of  converging 
cylindrical  and  spherical  magnetohydrodynamic  shocks  are  independent 


of  the  applied  field  and  are  given  by  the  usual  gas  dynamic  theory. 
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6*  Evaluation  of  F(a- ,  m^) 


Since  K  is  given  in  terms  of  v ,  the  integrated  form  of  the  shock 
Strength-area  relation  in  (5.  3)  is  that  chosen  for  numerical  evaluation. 

F  defines  a  one  parameter  family  of  functions  with  m^  as  parameter  and 
tables  of  this  function  for  various  values  of  m^  are  presented  in  the 
Appendix.  The  tables  refer  to  the  evaluation  of 


log  F((r  ,  m  )  =  / 

®  ^  1. 01 


0^-1  +  m^  y[  30-1  -  60 X  +  3(1+  e)x^  -  2x^]/2 


(e-x)(x-l)[  e+l  +  ym^  (x-1)V2]K2(x,  m^) 


dx 


The  tables  may  be  used  in  several  ways,  e.  g. ,  suppose  there  is  a 
section  of  variable  area  connecting  two  channels  of  constant  but  unequal  area, 
say  Aj  and  A^ ,  then  for  a  given  m^ , 

F(<r  ,  m  )  A 

If  Aj,  A^  and  <r^  are  known,  cr^  may  be  obtained  by  interpolating  inversely 
from  the  tables. 

After  passage  through  the  transition  section,  the  shock  travels  with  this 
altered  strength,  <r^,  and  ultimately  becomes  uniform.  Conversely,  given  o-^^, 
0-^  ,  m^  and  A^,  the  tables  may  be  utilized  to  determine  what  A^  would  give 
the  desired  o-^  .  Equation  (  6. 1)  could  also  be  considered  as  one  for  m^  with 
(Tj,  <r^,  Aj  and  A^  prescribed. 

Converging  cylindrical  and  spherical  shocks  may  be  treated  by  choosing 

2 

channels  with  areas  proportional  to  D  and  D  ,  respectively,  where  D  is 


#310 


-25 


the  distance  from  the  center. 

The  author  wishes  to  thank  Mr.  J.  Al-Abdulla  for  carrying  out  the 
computations  employed  in  this  paper. 
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The  tables  refer  to  the  evaluation  of  Z(o- ,  m^)  * 

2  2$  * 

<r  (0  ‘1) +m  y[ 30-1- 6ex +  3(1  + 0)x^-2x^]/2 
-  log^  F«r  ,  m^)  =  /  - i - - - !_ 

l.Ol[{0-x){x.l)[0  +  l+ Ymj{x-l)V2]K2((r,  m^) 

The  numbers  in  parentheses  refer  to  the  power  of  ten  by  which  each  entry  is 
to  be  multiplied. 


cr 

Z(<r.  0) 

Z(<r,  1) 

1.2 

6.65000 

(0) 

6.66275  (0) 

1.4 

1.36500 

(1) 

1.37038  (1) 

1.6 

2.  06500 

(1) 

2.07765  (1) 

1.8 

2.  76500 

(1) 

2.78834  (1) 

2,  0 

3.  46500 

(1) 

3.50263  (  1) 

2.2 

4.  16500 

(1) 

4.22054  (  1) 

2,  4 

4.86500 

(1) 

4.94203  (  1) 

2.  6 

5.56500 

(1) 

5.  66694  (  1) 

2.  8 

6. 26500 

(1) 

6.  39504  (  1) 

3.  0 

6. 96500 

(1) 

7.12610  (1) 

3.2 

7. 66500 

(1) 

7.85988  (  1) 

3.  4 

8.  36500 

(1) 

8.59614  (  1) 

3.  6 

9. 06500 

(1) 

9.  33473  (  1) 

3.  8 

9. 76500 

(1) 

1.  00755  (2) 

4.  0 

1. 04650 

(2) 

1.  08185  (  2) 

4.2 

1. 11650 

(2) 

1.  15638  (  2) 

4.  4 

1. 18650 

(2) 

1.23116  (2) 

4.  6 

1. 25650 

(2) 

1.30623  (  2) 

4.  8 

1. 32650 

(2) 

1.  38167  (2) 

5.  0 

1. 39650 

(2) 

1.45761  (  2) 

5.2 

1. 46500 

(2) 

1.  53426  (2) 

5.  4 

1. 53650 

(2) 

1.61206  (2) 

5.  6 

1. 60650 

(2) 

1.69194  (2) 

5.  8 

1. 67650 

(2) 

1.77700  (2) 
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Z(<r.  2) 


Z(<r.  3) 


6.70180  (0) 
1.38736  (1) 
2.  11842  (  1) 
2.  86366  (  1) 
3.62160  (1) 

4.38964  (1) 
5.  16488  (  1) 
5.94468  (  1) 
6.72696  (  1) 

7.  51020  (  1) 

8.  29343  (  1) 
9.07609  (  1) 
9.85792  (  1) 
1.  06390  (2) 
1.14194  (2) 

1.  22000  (2) 
1.  29803  (2) 
1.  37623  (2) 
1.45467  (2) 
1.  53355  (2) 

1.61318  (2) 
1.69410  (2) 
1.  77755  (  2) 
1.86748  (2) 


6.76868  (0) 
1.41736  (1) 
2.  18986  (1) 

2.  98974  (  1) 

3.  80758  (1) 

4.63410  (1) 
5.46244  (1) 
6.28849  (1) 
7.11011  (1) 
7.92647  (1) 

8.  73751  (  1) 

9.  54358  (1) 
1.03452  (2) 
1.  11432  (2) 
1.  19383  (2) 

1.  27314  (2) 
1.  35234  (2) 
1.43155  (2) 
1.  51093  (2) 
1.  59070  (2) 

1.67117  (2) 
1.  75300  (2) 
1.  83738  (2) 
1.  92857  (2) 
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Figure  1 

A  graph  of  the  parameter  IK  v's.  o-  for  m,  =  0 , 

U  1 

the  ordinary  gas  dynamic  case,  and  for  m^  =  2  and  6  . 


A  graph  of  the  parameter  vs.  o-  for  m^  =  100. 


I  2  3  4  5  6  g- 


Figure  4 

A  graph  of  the  parameter  K  vs.  it  for  m,  =  1000, 

Near  the  minimum  point  of  the  graph,  this  curve  crosses  that  for  m  =  100 , 
though  it  lies  beneath  the  latter  for  all  other  (t  .  This  indicates  a  possible 
reversal  of  effects  for  very  large  applied  fields,  i.  e.  ,  an  Increase  in  m 
may  lead  to  an  increase  in  K  .  More  extensive  numerical  investigation's 
are  thus  indicated. 
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Figure  5 


A  graph  of  the  parameter  6  vs.  o* . 


Since  there  is  a  small  total  variation  with 


m^ ,  only  the  case 


m. 


is  given. 


-34 


#310 


REFERENCES 

1.  Chester,  W.  "The  Quasi-Cylindrlcal  Shock  Tube",  Phil.  Mag.  (7), 

45,  1293,  1954. 

2.  Chisnell,  R.  "The  Motion  of  a  Shock  Wave  in  a  Channel,  With  Applications 

To  Cylindrical  and  Spherical  Shock  Waves",  J.  Fluid  Mech. ,  286,  1957. 

3.  Courant,  R. ,  and  Friedrichs,  K.  "Supersonic  Flow  and  Shock  Waves", 
Interscience  Publ. ,  New  York,  1948. 

4.  Friedrichs,  K.  "Nonlinear  Wave  Motion  in  Magnetohydrodynamics",  Los 
Alamos  Report  2105,  1957. 

5.  Germain,  P. ,  and  Gundersen,  R.  "Sur  les  ^coulements  unidimensionnels 
d*un  fluide  parfait  S  entropie  faiblement  variable",  C.  R.  Acad.  Sci. , 

Paris,  241,  925,  1955. 

6.  Gundersen,  R.  Thesis,  Brown  University,  Mar. ,  1956. 

7.  Gundersen,  R,  "The  Flow  of  a  Compressible  Fluid  With  Weak  Entropy 
Changes",  J.  Fluid  Mech. ,  3-,  553,  1958. 

8.  Gundersen,  R.  "A  Note  on  Shock  Flow  in  a  Channel",  J.  Fluid  Mech. ,  4, 
501,  1958. 

9.  Gundersen,  R.  "A  Perturbation  Analysis  of  Shock  Flow  in  a  Nozzle", 

J.  Aerospace  Sci. ,  26,  763,  1959. 

10.  Gundersen,  R.  "The  Piston- Driven  Shock",  J.  Aerospace  Sci. ,  27,  467, 

1960. 

11.  Gundersen,  R.  "The  Propagation  of  a  Plane  Shock  Wave  Into  a  Moving 
Fluid",  J.  Aerospace  Sci. ,  28,  584,  1961. 


?310  _35. 

12.  Gundersen,  R.  "Shock  Propagation  in  Time-Dependent  Ducts",  J.  Aerospace 
Sci. ,  ^  748,  1961. 

1  3.  Gundersen,  R.  "Magnetohydrodynamic  Shock  Propagation  in  Non-Uniform 
Ducts",  MRC-TS  287,  1961. 

14.  Gundersen,  R.  "Cylindrical  and  Spherical  Shock  Waves  in  Monatomic 
Conducting  Fluids",  MRC-TS  299,  1962. 

15.  Mirels,  H.  "Source  Distribution  for  Unsteady  One- Dimensional  Flows 
With  Small  Mass,  Momentum,  and  Heat  Addition  and  Small  Area 
Variation",  NASA  Memo  5-4-59E,  1959. 

16.  Rosciszewski,  J.  "Calculations  of  the  Motion  of  Non-Uniform  Shock 
Waves",  J.  Fluid  Mech. .  8,  337,  i960. 

17.  V\/hitham,  G.  "On  the  Propagation  of  Shock  Waves  Through  Regions  of 
Non-Uniform  Area  or  Flow",  J.  Fluid  Mech. ,  4,  337,  1958. 


